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ABSTRACT 

Current cosmological analyses which use Type la supernova (SN la) observations combine SN sam¬ 
ples to expand the redshift range beyond that of a single sample and increase the overall sample 
size. The inhomogeneous photometric calibration between different SN samples is one of the largest 
systematic uncertainties of the cosmological parameter estimation. To place these different samples 
on a single system, analyses currently use observations of a small sample of very bright flux standards 
on the HST system. We propose a complementary method, called ‘Supercal’, in which we use mea¬ 
surements of secondary standards in each system, compare these to measurements of the same stars 
in the Pan-STARRSI (PSI) system, and determine offsets for each system relative to PSl, placing all 
SN observations on a single, consistent photometric system. PSl has observed Stt of the sky and has 
a relative calibration of better than 5 mmag (for ^ 15 < griz < 21 mag), making it an ideal reference 
system. We use this process to recalibrate optical observations taken by the following SN samples: 

PSl, SNLS, SDSS, CSP, and CfAl-4. We measure discrepancies on average of 10 mmag, but up to 35 
mmag, in various optical passbands. We find that correcting for these differences changes recovered 
values for the dark energy equation-of-state parameter, w, by on average 2.6%. This change is roughly 
half the size of current statistical constraints on w. The size of this effect strongly depends on the 
error in the B — V calibration of the low-z surveys. The Supercal method will allow future analyses 
to tie past samples to the best calibrated sample. 


1. INTRODUCTION 

Since the initial discovery of cosmic accelerat ing ex¬ 
pansion (Riess et al.|1998 Perlmutter et al.|1999), many 
samples of Type la supernovae (SIN la) have been ac¬ 
quired to better constrain the dark energy equation-of- 
state parameter, w. As the systematic uncertainties in 
joint samples are nearly equal to the statistical uncer¬ 
tainties, increasing effort must be expended on reducing 
systematic uncertainties so that the total uncertainties 
do not soon hit a systematic floor. Of all the systematic 
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uncertainties, recent analyses (e.g., Betoule et al. 2014 
Scolnic et al.[2014a| [hereafter S14]) found that those re 


fated to photometric calibration make up > 70% of the 
total systematic uncertainty and pose the most immedi¬ 
ate challenge. 

Most SN analyses that attempt to constrain w combine 
publicly available SN samples to improve statistics and 
cover a wider redshift range. To reproduce the most re¬ 


cent c osmological results (Betoule et al.|2014 


20141, one must combine SIN from >10 independently 
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calibrated photometric systems. The calibration of each 
system is performed by multiple groups using sometimes 
significantly different methodology. To date, no cosmol¬ 
ogy analysis with SN la finds a solution so that the cal¬ 
ibration of all the various systems is consistent. With¬ 
out doing so, a SN analysis likely underestimates the 
systematic uncertaint ies of the cosmological parameters. 


Betoule et al. (2014| showed that there are differences 
of 5% between average distances of a single sample rel¬ 
ative to the expected distances from the ACDM model. 
It is now imperative to determine whether the small de¬ 
viations are due to noise, calibration uncertainties, or 
deviations from the ACDM model. 

In this analysis, we take advantage of the uniform cal¬ 
ibration of_the_^PSl_smvey_with <1% precision and ac¬ 
curacy (Schlafly et al. 2012) in order to measure and 
improve the consistency ot catalog photometry between 
different surveys. Here, we use publicly available data 
from HST, SDSS, SNLS, CSP, and CfAl-4. In O we 
discuss how each sample is currently calibrated. We in¬ 
troduce an analysis of the cross-calibration between mul¬ 
tiple systems in In Fj we discuss the magnitude of 
the discrepancies and the implications of correcting for 
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these, including effects on SNIa distances and recovered 
cosmology. We also quantify the dominant uncertainties 
in this approach. Our discussion and conclusions are in 
^and ^ 

2. CURRENT CALIBRATION 

In recent cosmological analyses, the calibration of each 
system may be compared d i rectly or indirectly to the AB 
system (Oke & Gunn'1983 Fukiigita et al.|1996 |. In the 
AB system, a monochromatic magnitude is dehned such 
that 


wab(j^) = -2.5 log 


10 




1 Jy 


8.90 mag, (1) 


where /i, is the flux per unit frequency from an object 
in Jy. Therefore, a magnitude 0 object should have the 
same counts as a source of fi, = 3631 Jy. 

We can define an AB broadband magnitude by the fol¬ 
lowing equation: 


ttiab = 2.5 X logio 


/ jhty) ^p{v)fu dv 
J{hv)~^p{v) 3631 Jy dv 


( 2 ) 


where p{v) is the filter response function and this equa¬ 
tion assumes that the detector is a photon-counting de¬ 
vice. 

Some of the surveys analyzed here calibrate their 
photometry to the AB system while others use Vega- 
calibrated systems. In either case, an AB offset can be 
given to convert the zeropoint of the calibration of one 
measured system to the true AB system. For a given 
system S and passband p, this may be expressed as 


S,p 

TTT-Ag = rn, 


S.p 

Sys 


aS,p 

^AB- 


( 3 ) 


where is the system magnitude, uiab 

magnitude and is the offset for a particular filter 

between the system magnitude and the AB magnitude. 

By Eq. 2, these offsets can be found explicitly if given a 
spectrum fi, defined to be on the AB system, a measured 
passband Pi, and the observed magnitude of the star in 
system S. For all systems below, the spectra u sed for 
this p rocess are taken of HST Calspec standards (Bohlin 
19961. These spectra are composite spectra from b'i'lb 
and NICMOS observations and have an uncertainty of 
^5 mmag for every 5000 A (Bohlin 2014) from 3000 A 
to 15000 A. 

There are two major components to the systematic un¬ 
certainties of each system’s calibration: how well obser¬ 
vations are tied to fundamental photometric standards 
and how well the photometric standards themselves are 
calibrated. The advantage of the calibration method for 
each survey discussed below is that it relies on obser¬ 
vations of known HST Calspec standards for which we 
have accurate spectra and can compare synthetic pho¬ 
tometry with observational photometry. The disadvan¬ 
tage of this approach is that the Calspec standards are 
sparse (^20 over the observable sky) and typically much 
brighter (r < 13 mag) than normal survey stars (r > 15 
mag). Furthermore, different surveys use different Cal¬ 
spec standards. 

Below, we briefly review the calibration of each of the 
photometric catalogs used in this analysis. We only an¬ 


alyze publicly available sources. A summary of the fil¬ 
ters, calibration standards, and systematic uncertainties 
of each system is given in Table 1. The systematic uncer¬ 
tainties for each system are described such that there is 
an uncertainty in the common zeropoint to observations 
of a given filter and an uncertainty in the mean wave¬ 
length of that filter. The transmission functions of all 
the filters from the various systems analyzed are shown 
in Fig. [2 

This section is separated into an explanation of the cal¬ 
ibration of the PSl photometric system, the systems of 
the other higher-z surveys and systems of the low-z sur¬ 
veys. The PSl, SNLS and SDSS systems share a common 
calibration path in that multiple Calspec standards are 
used to define the photometry on the AB system and 
nightly photometry can be tied directly to stellar cata¬ 
logs in each survey’s natural system. The low-z systems 
are partly tied to the Vega system and partly tied to the 
AB system, and only BD-f 17°4708 (hereafter BD+17) is 
used to tie magnitudes from each filter to the AB sys¬ 
tem. The AB offsets for all surveys are given in Table 
1. These offsets are the same ones applied when fitting 
light curves of the SN. 

For the low-z surveys, zeropoints of the nightly pho- 
tometry are determined e ither by transforming L andolt 
jLandolt fc Uomoto 20071 and Smith standards (Smith 
et al.||20U2p onto the respective natural systems or by 
transforming the nightly photometry onto the system of 
Landolt standards. The systematic uncertainties of the 
zeropoints given in Table 1 should therefore include: un¬ 
certainties in the AB magnitudes of the primary stan- 
dard(s) used, uncertainties in the measurements of the 
primary standard(s) by each survey, uncertainties in the 
transfer of zeropoints between the local standards to the 
primary standards, and systematic uncertainties in the 
measurements of the local standards. 

2.1. Calibration of the PSl full sky sample 

PSl - The PSl pho tometric calibration is presented 


Tonry et al. (2012b), hereafter T12. PSl is a 1.8 m 


telescope on Haleakala with a field of view of 7 square 
degrees. The observation strategy has two large parts: a 
Stt survey across the entire observable sky and a Medium 
Deep survey for 10 fields of 7 square degrees each. Both 
surveys observe in gpirpiipizpiypi (see Figure 1 for 
gpiTpiipiZpi fi lters used for this an alysis) . T12, based 
on wor k from Stubbs et al. (2010) and Tonry et al. 
( |2012aD , uses an innovative laser diode system to accu¬ 
ratelyand precisely determine the filter bandpass edges 
(a\ < 7 A) and throughput curves. The flux calibra¬ 
tion of PSl measurements relies o n an iterative process 
that includes work from T12 and Schlafly et al. (2012) 
and is augmented by S14. T12 analyzes observations of 
7 HST Calspec standards with PSl and compares the 
observed magnitudes of these standards to the predicted 
magnitudes from synthetic photometry. T12 then finds 
the AB offsets (Eq. 3) so that the observed magnitudes 
best matches the synthetic photometry, given fixed con¬ 
straints from measurements on the bandpass edges and 
shapes. 

It is necessary that photometry from all observations 
from the Att survey and Medium Deep survey can be 
linked by a single ze ropoint for each filter. To do so, 
Schlafly et al. (2012) uses the ‘Ubercal’ process ( Pad-| 
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TABLE 1 

Previously reported calibration differences between systems 


System 

Filters 

Standards 

Stan. Observation 

ZP Err 
[mmag] 

Wave Err 
[nm] 

AB offsets 
[mmag] 

Ref. 

PSl 

griz 

7* 

PSl 

[12,12,12,12] 

[0.7,0.7,0.7,0.7] 

[-23,-33,-24,-24] 

T12,S14,S15 

SNLS 

griz 

3* 

SNLS 

[3,6,4,8] 

[0.3,3.7,3.1,0.6] 

[-7,-8,-13,4.7] 

B13 

SDSS 

ugriz 

3* 

SDSS PT 

[8,4,2,3,5] 

[0.6,0.6,0.6,0.6,0.6] 

[-,-28,-14,-27,-20] 

M8,D10,B13 

CfAl/2 

UBVRI 

BD-l-17 

Landolt 

[100,15,15,15,15] 

[2.5,1.2,1.2,2.5,2.5] 

[-,131,6,-168,-410] 

Land 

CSP 

ugri.BV 

BD-l-17 

Smith Landolt 

[23,9,8,7,8,8] 

[0.7,0.8,0.4,0.2,0.7,0.3] 

[-,10,-4,-13,102,6/-7/2] 

Sll 

CfAK 

ri.UBV 

BD-l-17 

Smith Landolt 

[25,7,31,11,7] 

[0.7,0.7,2.5,0.7,0.7] 

[-3,-9,-,80,6] 

H09a 

CfAS 

UBVRI 

BD-l-17 

Landolt 

[70,11,7,7,20] 

[2.5,0.7,0.7,0.7,0.7] 

[-,71,7,-179,-406] 

H09a 

CfA4 

ri.BV 

BD-l-17 

Smith Landolt 

[31,11,7,25,7] 

[2.5,0.7,0.7,0.7,0.7] 

[-3/1,-9,80/135,6] 

H09b 


Summary of various sytems used in this analysis. The columns are: Filters used for observations, Standards used to determine the 
absolute flux zeropoints, Standard Observations that details with what telescope/camera the standards were originally observed with, 
zeropoint error claimed by the survey, wavelength error of the filter bandpasses, AB offsets to transform from the system magnitudes 
to AB magnitudes, and primary reference. The standards used by SNLS and SDSS are G191B2B, GD153, and GD71. The standards 
used by PSl are given in Table 2. AB offsets for the low- 2 : systems are generally large as the system was not defined on the AB system. 
Multiple values for certain offsets indicate multiple periods of a survey where the filters changed. 
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Fig. 1. — The filter transmission functions of all the systems in 
this analysis. Comparisons are broken into griz and BVRI. For 
both CSP and CfA-Keplercam (CfAK), multiple filters are shown 
as in different periods of the survey, different filters were used. 


manabhan et al.|2008 ) that creates a relative calibration 
across the sky and directly ties the zeropoints of all pho¬ 
tometry to the catalogs from the fields that contain the 
observed Calspec standards analyzed in T12. Because 
PSl repeatedly observed the same regions ^12 times in 
each filter and there are overlapping regions between dif¬ 
ferent pointings, all relative zeropoints can be determined 
simultaneously and robustly. This procedure determines 
the system throughput, atmospheric transparency, and 
large-scale detector flat field. The solution from this 
process includes data from both the Medium Deep fields 
and the Stt sky over the full survey. From Ubercal, new 
photometric catalogs are created across the entire ob¬ 


servable sky with relative accuracy and precision better 
than 5 mmag. To improve on the initial zeropoints given 
by T12, S14 then iterates on this process by analyzing 
the entire sample of Calspec standards observed over the 
course of the Stt survey and redetermines the AB offsets 
for the PSl system. S14 also analyzes the variation of 
filter transmission functions across the focal plane and 
finds the variations across the focal plane with a radial 
dependence for stars with 0.4 < gpi — zpi < 1.5 to be 
up to 8 mmag (due to the design of the filters), though 
typically at a dispersion level of ^2 mmag. The filter 
throughput is measured at multiple radial positions, and 
brightnesses can be corrected based on radial position on 
the focal plane. 

For the present analysis, we repeat the process done 
in S14 to redetermine the AB offsets for the PSl sys¬ 
tem. We find that the majority of the observations of 
the Calspec standards from T12 placed the standards 
at the same position on the same chip for each obser¬ 
vation. This position was very close to the center of 
the focal plane where it has been noted tha t there is a 
stron g gradient in the behavior of the chip (Rest et al. 
20141. T he Ubercal solution for the large-scale detector 
hat field ( Schlafly et al.|20i2 ) at the location where these 
standards were observed varies by up to 20 mmag over 
less than a quarter of a CCD. Given this issue, observa¬ 
tions from T12 are not included in the present analysis. 
Also, only observations fainter than the saturation limit 
of [14.3, 14.4, 14.6, 14.1] in (/pirpiipizpi are included. 
Table 2 shows both the synthetic photometry of the six 
Calspec standards and the Ubercal photometry of these 
standards. As shown in Table 3, we find that corrections 
from S14 of 20-35 mmag in each filter are needed. The 
significant size of these corrections is partly due to the 
update of the HST Calspec spectra. Further information 
will be given in the public release of the PSl data. 

2.2. Calibration of intermediate and high-z surveys 

SDSS-II - The basis for the Sloan Digital Sky 
Survey-II (h ereafter referred to as ‘SDSS’) calibration is 


presented in Holtzman et al. (2008 ) and Doi et al. ( 2010 ), 
with important updates given m Betoule et al. (2013). 
The primary instrument o f the SPSS Supern ova Survey 


is the SDSS CCD camera (Gunn et al.||1998 ) which was 
moun ted on a dedicated 2.5 m telescope (jGunn et al. 
2006) at Apache Point Observatory (APO), blew iVlex- 
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TABLE 2 

PSl Observed and Synthetic Magnitudes of Calspec 

STANDARDS 


TABLE 3 

PSl Photometric System 


Star 

Filter 

Obs. Magnitude 

Syn. Magnitude 

SF1615 

gpi 

16.975 (0.007) 

16.996 

Snap-2 

5pi 

16.413 (0.008) 

16.447 

Snap-1 

5pi 

15.477 (0.007) 

15.495 

WD1657-I-343 

ffpi 

16.212 (0.007) 

16.224 

KF06T2 

ffpi 

14.391 (0.007) 

14.429 

lds749b 

ffpi 

14.562 (0.010) 

14.573 

C26202 

ffpi 

16.651 (0.008) 

16.676 

SF1615 

rpi 

16.523 (0.007) 

16.562 

Snap-2 

rpi 

16.012 (0.008) 

16.046 

Snap-1 

rpi 

15.857 (0.007) 

15.894 

WD1657-I-343 

rpi 

16.669 (0.007) 

16.693 

KF06T2 

rpi 

13.573 (0.010) 

13.606 

lds749b 

rpi 

14.765 (0.008) 

14.808 

C26202 

rpi 

16.347 (0.007) 

16.368 

SF1615 

ipi 

16.360 (0.006) 

16.385 

Snap-2 

ipi 

15.878 (0.007) 

15.904 

Snap-1 

ipi 

16.191 (0.007) 

16.202 

WD1657-I-343 

ipi 

17.053 (0.006) 

17.073 

lds749b 

ipi 

15.000 (0.010) 

15.039 

C26202 

ipi 

16.240 (0.008) 

16.263 

GD153 

2pi 

14.230 (0.007) 

14.263 

P177D 

2pi 

13.137 (0.008) 

13.154 

SF1615 

2pi 

16.285 (0.006) 

16.318 

Snap-2 

2P1 

15.846 (0.006) 

15.875 

Snap-1 

2pi 

16.393 (0.006) 

16.424 

WD1657-I-343 

2pi 

17.346 (0.007) 

17.360 

KF06T2 

2P1 

13.083 (0.010) 

13.084 

C26202 

2P1 

16.211 (0.006) 

16.245 

GD153 

ypi 

14.450 (0.007) 

14.472 

P177D 

ypi 

13.135 (0.007) 

13.135 

SF1615 

ypi 

16.269 (0.006) 

16.278 

Snap-2 

ypi 

15.837 (0.007) 

15.853 

Snap-1 

ypi 

16.581 (0.007) 

16.567 

WD1657-I-343 

ypi 

17.570 (0.006) 

17.579 

KF06T2 

ypi 

12.980 (0.007) 

12.991 

lds749b 

ypi 

15.380 (0.008) 

15.401 

G26202 

ypi 

16.225 (0.007) 

16.251 


The calspec standard stars with adequate PSl photometry are 
presented here. Both the observed and synthetic magnitudes of 
these standards are given, before an AB correction is applied to 
the natural PSl magnitudes (from Tonry et al. 2012). The 
uncertainties in these measurements are given in brackets. The 
synthetic spectra can be found on the Calspec websit^ we use 
version 005 in this analysis. 

‘‘http://www.stsci.edu/hst/observatory/crds/calspec.html 


ico. The survey observ ed in the five optical bands: ugriz 
(Fukugita et al. 1996 see Figure 1 for griz filters used 
for this analysis). I'he survey covers a 300 square-degree 
region (2.5° wide over 8 hours in right ascension). The 
photometry explored in this analysis is of the region cen¬ 
tered on the celestial equator referred to as ‘Stripe 82’. 
The meas urement o f SPS S effective passbands are de¬ 
scribed in |Doi et al.| (|2010|). 

The absolute flux calibration has been determined us¬ 
ing the SDSS Photometric Te lescope (PT) observ ations 
of Calspec solar analog stars ( Tucker et al. 2006 ) . Thi s 
process is described in detail in Holtzman et al. (2008). 


Betoule et al. (2013) updated the AB offsets to reflect 
recent revisions to the HST Calspec observations and 
SDSS observations. Small differences between the indi¬ 
vidual filters mounted on the different columns of the 
camera can be neglected due to the SDSS calibration 


Filter 

S14 

S15 

S14-S15 

S14-S15 


[mmag] 

[mmag] 

HST-Recal 

Improved Phot 

gpi 

-8± 12 

-20 ±8.0 

-7 

-13 

rpi 

-10 ± 12 

-33 ±8.0 

-9 

-24 

ipi 

-4± 12 

-24 ±8.0 

-9 

-15 

2pi 

-7± 12 

-28 ±8.0 

-9 

-12 


Corrections of the AB offsets from the original definition of the 
PSl calibration as given in in T12. The S15 include the updates 
to the most current HST Calspec magnitudes (version 005). The 
breakdown of the change in values due to the HST version update 
and from our own improved photometry is given in the last two 
columns respectively. The list of Calspec standards and both their 
synthetic and observed magnitudes used for this recalibration is 
given in Table 2. 


strategy (|Betoule et al.||2013|). 

SNLS - The latest calibration of the Supernova 


Legacy Survey (SNLS) is given in Betoule et al. 
SNLS uses the 3.6 m CFHT atop Mauna Kea. 


(2013). 

^NTS 


covers the four low extinction fields of the CFHT Legacy 
Survey Deep component (called D1 to D4). The fields are 
repeatedly imaged in the 4 optical bands gnA m 
(see Figure 1). The original im filter was broken in July 
2007 and replaced by a slightly different jm (denoted 
i2jvf) filter in October of the same year. As im mag¬ 
nitudes are given for all SN observations published, and 
not i2M, we exclude z2m from this analysis. Th e field-of- 
view of the MegaCam camera is 0.96x0.94 deg^ (Boulade 
et al.||20(53 ). Measurements of the SNLS bandpasses are 
presented in Regnault et al. (2009). 

The SNLS calibration is determined from multiple 
paths including observations of HST Calspec standards 
and Landolt stars. Spatial variations of the passband 
response result in variation in the brightness found for 
stars when they are observed at different positions on 
the MegaCam focal plane (R09). In the SNLS data re¬ 
lease, magnitudes of the stars are transformed as if they 
were observed at the center of the focal plane. 

2.3. Calibration of low-z surveys 


In recent SN analyses (Betoule et al. 2014 S14), the 
calibration of the absolute hux of the various low -2 sur¬ 
veys is tied to measurements of the primary standard 
star BD-l-17. For filter s ugri, the flux is calibrated to 
the Smith et al. (2002) magnitudes of BD-l-17. For fil¬ 
ters (7nF777~)~tnenujris calibrated t o the Landol t mag- 


nitudes of BD-l-17. Magnitudes from Smith et al. (2002) 
are expected to b e consistent with the AB sy stem at bet- 
ter than 4 mmag. Landolt & Uomoto (20071 showed that 
the Landolt magnitudes of various standard stars and the 
AB magnitudes are consistent to 6 mmag. 


Data from the Calan/Tololo survey (Hamuy et al. 
19931) are not included in this analysis as the published 


number of comparison stars is quite small and a large 
fraction of the stars are below the —30° declination limit 
of the PSl 37r survey. 

CSP - The basis for the CSP calibration is presented 


in Contreras et al. (2010). The CSP optical follow-up 
campaigns were carried out with the Direct CCD Cam¬ 
era attached to the Henrietta Swope 1 m telescope lo¬ 
cated at the Las Campanas Observatory (LCO). The 
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survey observed in ugriBV and the field-of-view of the 
observations is 8.7' x 8.7'. Definitive measurements of 
the CSP filter throughput curves were carried out at the 
telescope using a monochro mator and calibrated p hoto- 
diodes (Rheault et al. 2010; |Stritzinger et al.||20li|). 

CSP SN magnitudes are published in the native photo- 
metric system, defined by the S wope filter response func¬ 
tions of |Stritzinger et al.||2011 and the primary standard 
BD-l-17. The CSP local standard magnitudes are pub¬ 
lished in the standard system, and we convert these mag¬ 
nitudes onto the CSP native system using the transfor- 
mation equa tions provided. As discussed in |Stritzinger| 


et al. (2011), we use the three different filter transmis- 
sion functions for the CSP V band (shown in Fig. 1) for 
the three periods of the CSP survey in which a different 
filter was used (labeled ‘CSPb, ‘CSP 2 ’, ‘CSP 3 ’ for this 
analysis). 

CfA4 - The basis fo r the CfA4 calibration is presented 


Hicken et al. (2012). The CfA4 data were obtained on 


the 1.2 m telescope at the Fred Lawrence Whipple Ob¬ 
servatory (FLWO) using the single-chip, four-amplifier 
CCD KeplerCam4. Observations were acquired in a field 
of view of approximately 11.5' x 11.5'. Cramer et al. 
(in prep) measured the FLWO 1.2 m KeplerCam BVr'i' 
passbands using the mono chromatic illumination tech- 
nique initially described in Stubbs & Tonry (2006). No 
atmospheric component is included in the Keplercam fil¬ 
ter tra nsmission curves. As discussed in IHicken et al.l 
(2012), the 1.2 m primary mirror deteriorated during the 
course of the CfA4 so that different transmission func¬ 
tions for various filters were recognized for different parts 
of the survey (labeled ‘CfA4_l’ and ‘CfA4_2’ in this anal¬ 
ysis). 

CfA3 - The basis for the CfA3 calibration is pre¬ 
sented in Hicken et al. (|2009 ). The CfA3 sample was 
acquired on the F. L. Whipple Observatory 1.2 m tele¬ 
scope, mostly using two cameras, the 4Shooter camera 
and Keplercam. The field of view of the observations 
was approximately 11.5' x 11.5'. UBVRI filters were used 
on the 4Shooter (hereafter referred to as ‘CfAS’) while 
UBVr'i' filters were used on Keplercam (hereafter re¬ 
ferred to as ‘CfAK’). The 4Shooter BVRI passbands are 
given in Jha et al. (|2006 ). The Keplercam UBVr'i' are 
measured in Cramer et al. (in prep.), as discussed in the 
CfA4 subsection above. 

CfA2 - The basi s for the CfA2 calibration is presented 


Jha et al. (20061. The CfA2 sample was acquired on 


the 1'. L. Whipple Observatory 1.2 m telescope, with ei¬ 
ther the AndyCam CCD camera or the 4Shooter camera. 
Like for CfA3, the field of view was 11.5' x 11.5'. UBVRI 
filters were used on both cameras. 

CfAl - The basis for the CfAl calibration is presented 


Riess et al. 


acquired on t. 


(1999). Like, CfA2, the CfAl sample was 
le 1. 


j. Whipple Observatory 1.2 m tele¬ 
scope, which has a field of view was 11.5' x 11.5'. BVRI 
filters were used for all observations. To include this sam¬ 
ple in our analysis, we created our own catalogs of the 
stellar photometry presented by matching the catalogs 
given in the paper to the finding charts in the paper. 


3. CROSS CALIBRATION 

Here we present a new calibration path: to directly tie 
all photometry catalogs to the homogeneous PSl-Ubercal 
catalog. The PSl catalog covers Stt of the sky with ~5 
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Fig. 2. — Visual representation of the steps of the Supercal 
method. The top panel shows the synthetic and observed differ¬ 
ences in brightness of stars observed with the gsDSS ppiband. 
Observed stars are shown in black and synthetic stars are shown 
in red. A linear transformation is fit to both the synthetic and 
observed trend. In the middle-left, a histogram is shown of the 
residuals from subtracting the synthetic trend from the observed 
stellar photometry. The histograms are normalized to unity. In the 
middle-right, the residuals as a function of magnitude are given. In 
the bottom panel, the dependency of the residuals on sky position 
is shown (bins given as yellow triangles). 


mmag relative calibration, and depth to ^22 mag. In 
our new method, we compare observed color differences 
of stars in common to surveys with the synthetic trans¬ 
formations between the systems. We then correct for 
these differences so that all samples are calibrated in a 
consistent manner. The advantage of this new method 
is that there are much better statistics as we can use 
many of the stars in the catalogs that are directly used 
to determine zeropoints of the SN photometry. 

For a given star in which photometry is measured by 
two surveys in the AB system, the expected difference in 
magnitude can be expressed as: 

Osi - Os2 = Ysi - Yg2 + A5 i _52 (4) 

p p' y p' p p' 

where Ogi and O52 are the observed magnitudes in pass- 

P p' 

band p and p' in the two systems, Ygi , Yg 2 are the syn- 

P p' 

thetic magnitudes (Eq. 2,3) of the systems and A5i_g2^ 

is the systematic discrepancy between how the calibra¬ 
tion systems are defined. Typically, A5i_52^ is assumed 

to be 0; the main goal of this analysis is to test this 
assumption. A5i_52^ 7^ 0 if there are errors in the mea¬ 
surements of the bandpass definitions, errors in the pho- 
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tometry of the star compared, errors in the synthetic 
spectrum of the star itself or errors in the AB offsets (Eq. 
3) given for each system. Given a number of standards 
with well-calibrated spectra, from Eq. 2 both A5i_52^ 

and wavelength shifts to p{iy) and p'(y) can be found by 
comparing the observed magnitudes of stars from two 
surveys. 

Our new approach consists of finding the differences of 
observed magnitudes of the same stars observed by two 
surveys, and comparing these differences with the ex¬ 
pected differences from known passband definitions and 
a stellar library of synthetic spectra. To do so, we fit a 
line to both the observed differences Osi — Og 2 of the 

P p' 

two systems with color and the expected synthetic dif¬ 
ferences hsi — Yg 2 of the two systems with color: 

P p' 

Osi - 0^2^ = Q!o(Osi - dgi^ - Co) -I- Po (5) 
= ctriYsi - Ygi^ - Co) -I- /?f (6) 

where O and Y describe the vector of observed and syn¬ 
thetic magnitudes of all overlapping stars between two 
surveys (input values), a and (3 are the linear compo¬ 
nent of the fit to the differences between the magnitudes 
for the observed and synthetic sequences (fitted values) 
and Co is a reference color that is chosen. The linear fit 
is performed versus a transformation color (O51 — 051^ 

or Y51 — Ygi^) as we expect differences in magnitude 
between surveys with similar passbands to depend on 
color. If the passband shapes and edges are accurate 
but there is still a discrepancy in the zeropoints, then 
Po — Py = A5 i_ 52^ at the reference color, co. If the 

passband shapes and edges of either system are incor¬ 
rect, then we should expect that ao — cay 0. It is 
possible that ao = ay or Po = Py if systems are con¬ 
sistently wrong, though this scenario is unlikely. 

To compare photometry from two different surveys, we 
follow a seven-step process in order to determine a single 
mean offset between the catalogs. 


magnitudes to determine (Ygi — Yg 2 ). Similarly, 

P p' 

find the color (Ygi — Ygi^). 

4. For the observed sequence of stars, adopt a mag¬ 
nitude cut to reduce the Malmquist bias from the 
PSl faint stars included in the sample. No magni¬ 
tude cut is determined from the other samples as 
all catalog stars are used in the external analyses 
to determine the calibrate the SN photometry. We 
also adopt a bright-end cut for PSl magnitudes of 
[14.8,14.9.15.1,14.6] in gp 1 rpiipizpibecause of con¬ 
cerns about linearity brig hter than these magni¬ 


tudes (Schlafly et al. 


2012 


TT 


5. For the observed sequence of stars, correct the stel¬ 
lar magnitudes for Milky Way reddening using the 
know n positions of the stars and e xtinction values 
from Schlafly & Finkbeiner (2011) with a 2D sky 
map. 'I'he extinction values are specific to each 
system. 


6 . Choose a specific color range (of e.g., g — i) oi the 
catalog stars used in the analysis to reflect the same 
colors of the stars in the synthetic library. 


7. Fit ao, Po and ay, Py to the observed and syn¬ 
thetic sequences respectively. All photometric er¬ 
rors are propagated and we perform iterative 3 — cr 
clipping to determine best fit values. 

An illustration of these steps is shown in Fig. The 
offset A5 i _52 is determined from comparing the differ- 

PP' 

ences in /3 at a reference color cq of the fitted lines to the 
synthetic and observed differences. The statistical errors 
in either case includes both the error in the measurement 
of the linear ht to the observed and synthetic distribu¬ 
tions. We check the dependence of the offset on both 
brightness and R.A. but only make a cut on brightness 
to prevent a Malmquist bias. We do not correct for a 
spatial bias any survey within this analysis, though this 
can be done in future analyses. 


4. RESULTS 


1. Match the astrometric positions of stars observed 
by two surveys to < 1 arcsec. This ensures that 
stars will not be mismatched. To avoid potential 
errors related to blending, only isolated stars, those 
with no other star with m < 22 mag within a 15 
arcsec radius, are included in the sample. 

2. For a given band (e.g., p), or bands from differ¬ 
ent systems that are near each other in wavelength 
space, subtract the observed magnitudes from the 
two matched catalogs {Ogi —Og 2 ). Also determine 

P p' 

the transformation color to use {Ogi —Ogi^ where 
p and p2 are two passbands chosen to provide the 
strongest and most linear leverage on the color tilt 
of the magnitude differences. 

3. For the chosen bands, integrate the spectral library 
through the passbands of both systems to deter¬ 
mine the synthetic magnitudes in those passbands. 
For these library spectra, subtract the synthetic 


4.1. Discrepancies Between Surveys 

We follow the procedure explained in the previous sec¬ 
tion to determine discrepancies between PSl and each 
available catalog. For each comparison, we ht the differ¬ 
ences between observations from two different systems 
but with similar passbands, and attempt to remove the 
dependence on color. Since the main purpose of this anal¬ 
ysis is to more consistently tie all of the systems to the 
HST Calspec system, we use the HST Calspec library 
as our spectral library to determine the color transfor¬ 
mation between the two systems. All Calspec standards 
used in this analysis have a relatively small color range 
of 0.35 < g — i < 0.55 mag. Since we do not extrapolate 
beyond this color range, the number of stars used in the 
comparisons is limited. This is further discussed in §4.2[ 
For the higher-z surveys, we choose the transformation 
color for these comparisons to be gpi —ipi as we hnd this 

Saturation is observed in PSl magnitudes at 
[14.3,14.4.14.6,14.1] in gpirpiipizpi, and we conservatively 
add a half magnitude here. 
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choice results in the smallest scatter in residuals. We find 
other choices for the transformation color produce sim¬ 
ilar results, but with larger uncertainties. Additionally, 
using gpi — zpi rather than a color from another survey’s 
observations allows us to be consistent in our compar¬ 
isons (see Figure 2). However, because of issues at blue 
wavelengths including the Balmer jump, in order to com¬ 
pare gpito the B band, we use the transformation color 
B — i. In this case, we must perform a two-dimensional 
minimization as we find discrepancies in B — gpi ver¬ 
sus B — ipi, where the B measurement comes from the 
comparison sam ple (see discussion of systematic uncer¬ 
tainties in j |4.2[ ). 

We present the results from the offsets (A s^-s ^,) found 

after removing the color transformation (see Fig. 2) in 
Fig. 3. Here, Agi_g 2 ^ = j3o ~ Py and we analyze the 


very small color range of HST Calspec standards such 
that ao — oiY is insignificant. Offsets with respect to 
gpiTpiipiZpi are shown for each filter of each system. 
The errors shown include systematic uncertainties which 
are discussed in the next section. We separate SNLS 
into its 4 deep fields and SDSS-H into the two parts of 
Stripe82 that overlap with PSl Medium Deep fields and 
one part for the entire Stripe82. The largest deviations 
are seen in the comparison between g and B bands. Scat¬ 
ter of these offsets is around 2 - 3%. The offsets rela¬ 
tive to the rpi and ipi bands are generally < 10 mmag 
though there appears to be a systematic offset between 
PSl and other surveys in the r band. The offsets found 
relative to z band are slightly larger; in particular, there 
is a ^15 mmag offset between PSl and SDSS. The errors 
shown in Fig. 3 account for both the uncertainty from the 
observed magnitudes of stars as well as the uncertainty 
in the transformation from the synthetic magnitudes of 
the library standards. The uncertainties are largest for 
the comparison between gpiand the low-z systems’ B 
because there is the most scatter in the transformation 
between these two filters and the number of stars used 
for these comparisons is typically not high. 

For SNLS, we are able to perform these comparisons 
for each of their deep fields independently and find scat¬ 
ter < 5 mmag which shows both the precision of the 
PSl Ubercal and of the SNLS calibration. Similarly, we 
can verify the relative calibration of the PSl Ubercal in 
the comparisons with SDSS for the two Medium Deep 
fields that overlap with Stripe82. The largest difference 
between the two Medium Deep fields when comparing 
against SDSS is 4 mmag. Overall the differences seen 
here between the SDSS and SNLS calibration is ~ 10 
mmag in each band, depending on the field. This is 
roughly within the erro rs given for the join t analysis be¬ 
tween the two surveys (|Betoule et al.||2013|). 

In the Table 4, we present the offsets of /io — Py using 
the HST Calspec library at cq = {g — i)o = 0.45 mag, and 
both Po — Py and ao — av when using the NGSL library 
and a large enough color range to measure the slopes of 
the color transformations. We also give the number of 
stars common to PSl and each survey for each compari¬ 
son. As discussed below, we cannot accurately quantify 
the systematic uncertainties of the measurements of the 
slope difference ao — aY so the information given here 
is solely for future study. However, we note the large 
magnitude of the differences in slopes for some of these 


comparisons. For better understanding, we convert the 
difference in slope to a nominal difference in the mean 
effective wavelength for the filter used in the given sur¬ 
vey. Both the dependence of the slope on mean effective 
wavelength as well as the corresponding shift necessary 
to bring the calibration into agreement with PSl is given. 
Many of these values are larger than the quoted system¬ 
atic uncertainties in Table 1. 

The largest discrepancies seen in the B band of the 
low-z systems are most likely due to the use of BD-l-17 
as the primary calibration standard. |Bohlin &: Landolt] 
(20151 show that the luminosity of this standard has var¬ 
ied over the last two decades by ~4% and using it to an¬ 
chor the calibration will result in additional systematic 
biases. If we instead recalibrate the low-z surveys based 
on the Calspec standyds P 177D or P330E, as suggested 
by|Bohlin & Landolt (20151, we find that the change in B 
magnitudes is negligible, but the change in VRI magni¬ 
tudes is —30 mmag. The net result would be a consistent 
30 mmag offset in every filter when comparing the low- 
z BVRI to the higher-z systems. Therefore it appears 
likely that there is indeed an error in the low-z systems 
of A{B — U) = 30 mmag and possibly an additional gray 
offset of a similar magnitude. 

One can adjust the PSl zeropoints such that the 
PSl calibration is aligned to any other system. Do¬ 
ing this, we can then measure the discrepancies be¬ 
tween the recalibrated PSl system and all other sys¬ 
tems. For example, correcting the PSl system to align 
with SNLS requires applying offsets of —12,—6,-1-1,-1-7 
mmag from gpirpiipiZpi, respectively. We can then de¬ 
termine the discrepancies between all other systems and 
the PSl-l-SNLS system. 

We release onlincF^ all of the data used to make each 
comparison, and all iterations of versions of Figure 2 for 
all of the comparisons done in this analysis. We also 
include online all of the transmission functions for each 
of the systems used in this analysis, as well as all of the 
zeropoints for each filter. 


4.2. Systematic Uncertainties 

The largest systematic uncertainties in this approach 
are: the consistency of the PSl Ubercal, accuracy of the 
spectral library, reddening of the stars from dust and 
PSl linearity. A smaller systematic uncertainty is the 
variation of the filter functions with radial position on the 
focal plane, though for both PSl and SNLS, knowledge of 
this variation is used to correct the magnitudes based on 
their radial position. For all other systems, this variation 
is expected to be negligible. A summary of the dominant 
systematic uncertainties is given in Table 5. 

The systematic uncertainties from the PSl Ubercal are 
given in Finkbeiner et al. (2015, submitted). By compar¬ 
ing the relative calibration of SDSS and PSl, Finkbeiner 
et al. 2015 finds that for a given sky position, there 
are systematic uncertainties less than 9,7,7,8 mmag in 
griz. These uncertainties include both the systematics 
from SDSS and PSl, not only PSl. A fair upper limit 
for the systematic uncertainty due to PSl is roughly half 
these uncertainties: 5,4,4,4 mmag. These uncertainties 
are further reduced for the PSl Medium Deep fields due 
to the large number of observations of these fields and 


www.kicp.uchicago.edu/~dscolnic/supercal/ 
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Fig. 3. — The deviations from agreement with PSl calibration for all surveys and filters. The surveys shown are CfAl, CfA2, CfAS, 
CfAK, CSP, CfA4, SNLS, SDSS. The system and filter are listed for each comparison.. The uncertainties shown here include both the 
statistical and systematic uncertainties from the Supercal process. The offsets shown here should be added to each system’s magnitudes to 
agree with PSl. 
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TABLE 4 

Calibration discrepancies with PSl 


Survey 

Filtl 

[PSl] 

Filt2 

[O] 

NStar 

HST Off. {(S) 
(mmag) 

NGSL Off.* (/3) 
(mmag) 

Slopeobs (o) 
(mmag/mag) 

Slopes.y„ {a) 
(mmag/mag) 

Slope Diff. 
(mmag/mag) 

da/dX 

/ mmag \ 

AA 

(nm) 

(mag X nm) ^ 

CfAl 

9 

B 

31 


-15.5 ± 18.6 

-302.2 ± 54.2 

-334.3 ±9.1 

32.1 ± 54.9 

9.4 

3.4 ± 5.8 

CfA2 

9 

B 

146 

-1.3 ± 13.5 

8.4 ± 13.4 

-290.4 ± 9.7 

-334.3 ±9.1 

43.9 ± 13.3 

9.4 

4.7 ± 1.4 

CfAK 

9 

B 

731 

27.2 ±5.3 

30.9 ±5.3 

-317.6 ± 8.0 

-315.3 ± 12.2 

-2.3 ± 14.6 

8.7 

-0.3 ± 1.7 

CfAS 

9 

B 

396 

30.8 ±5.3 

34.5 ±5.3 

-293.3 ± 9.0 

-314.3 ± 12.2 

21.0 ± 15.1 

9.4 

2.2 ± 1.6 

CfA4i 

9 

B 

855 

19.8 ±5.7 

20.1 ± 5.7 

-349.7 ± 10.0 

-307.4 ±23.5 

-42.3 ± 25.5 

8.7 

-4.9 ±2.9 

CfA42 

9 

B 

912 

-6.4 ± 6.4 

-5.2 ± 6.4 

-315.1 ± 9.3 

-247.9 ± 14.1 

-67.2 ± 16.9 

9.3 

-7.2 ± 1.8 

CSPl 

9 

B 

251 

4.6 ±4.7 

1.4 ±4.7 

-321.7 ± 12.5 

-299.2 ±20.7 

-22.5 ± 24.2 

12.8 

-1.8 ± 1.9 

CSPl 

9 

9 

362 

-12.5 ±3.9 

-11.9 ± 3.8 

-91.8 ± 10.8 

-64.5 ±3.2 

-27.4 ± 11.2 

5.8 

-4.7 ± 1.9 

SNLSO 

9 

9 

200 

-11.5 ±4.1 

-17.9 ±4.1 

8.2 ±5.3 

18.0 ± 1.3 

-9.8 ± 5.4 

8.0 

-1.2 ±0.7 

SNLSl 

9 

9 

272 

-11.4 ±3.7 

-16.8 ± 3.7 

18.0 ±6.4 

18.0 ± 1.3 

-0.1 ± 6.5 

8.0 

-0.0 ±0.8 

SNLSl 

9 

9 

264 

-10.4 ±6.5 

-14.7 ± 6.5 

14.4 ±6.4 

18.0 ± 1.3 

-3.6 ±6.5 

8.0 

-0.4 ± 0.8 

SNLS3 

9 

9 

568 

-11.5 ±6.4 

-19.5 ± 6.4 

12.5 ±5.4 

18.0 ± 1.3 

-5.5 ± 5.5 

8.0 

-0.7 ± 0.7 

SDSSO 

9 

9 

3300 

-0.8 ±4.4 

-0.6 ±4.4 

-115.8 ± 1.1 

-121.3 ±4.6 

5.5 ±4.7 

8.7 

0.6 ±0.5 

SDSSl 

9 

9 

1525 

-0.5 ±4.4 

-1.4 ±4.4 

-117.5 ± 1.5 

-121.3 ±4.6 

3.8 ±4.9 

8.7 

0.4 ± 0.6 

CfAl 

r 

V 

22 


-24.4 ± 10.3 

-325.9 ±46.1 

-314.8 ±3.8 

-11.1 ± 46.3 

5.4 

-2.1 ± 8.6 

CfA2 

r 

V 

120 

-6.6 ±6.6 

-11.4± 6.4 

-364.9 ± 14.8 

-314.8 ±3.8 

-50.1 ± 15.3 

5.4 

-9.3 ±2.8 

CfAK 

r 

V 

634 

-1.1 ± 3.6 

-2.6 ±3.6 

-357.4 ± 10.8 

-348.8 ±5.9 

-8.6 ± 12.3 

5.1 

-1.7± 2.4 

CfAS 

r 

V 

339 

2.1 ± 3.8 

1.8 ±3.8 

-322.3 ± 16.5 

-359.7 ± 6.3 

37.4 ± 17.7 

5.4 

6.9 ±3.3 

CSPl 

r 

V 

309 

-8.4 ± 3.6 

-12.6± 3.7 

-393.1 ± 11.5 

-362.0 ±6.6 

-31.1 ± 13.3 

5.3 

-5.9 ±2.5 

CSPa 

r 

V 

309 

-2.2 ±3.7 

-3.4 ±3.8 

-376.3 ± 11.4 

-376.8 ±6.6 

0.5 ± 13.2 

5.2 

0.1 ± 2.6 

CSPs 

r 

V 

309 

-6.2 ±3.9 

-9.2 ±4.0 

-395.1 ± 11.6 

-380.0 ±7.2 

-15.1 ± 13.6 

5.2 

-2.9 ±2.6 

CfA4i 

r 

V 

1038 

-5.4 ± 3.6 

-4.5 ± 3.6 

-355.8 ± 9.6 

-348.8 ±5.9 

-7.0± 11.2 

5.1 

-1.4 ± 2.2 

CfA42 

r 

V 

1038 

-5.9 ±3.7 

-5.1 ± 3.7 

-355.8 ± 9.6 

-346.4 ± 5.9 

-9.4 ± 11.2 

5.1 

-1.9 ±2.2 

CfAl 

r 

r 

19 


-4.1 ± 13.0 

177.7 ±46.4 

76.4 ± 1.9 

101.3 ±46.4 

2.7 

36.9 ± 16.9 

CfA2 

r 

r 

119 

5.2 ± 5.1 

3.6 ±4.8 

13.8 ±20.7 

76.4 ± 1.9 

-62.6 ± 20.8 

2.7 

-22.8 ± 7.6 

CfAS 

r 

r 

340 

13.5 ±3.9 

14.8 ±3.9 

93.3 ± 16.0 

79.7 ± 2.3 

13.6 ± 16.2 

2.7 

5.0 ± 5.9 

CfAK 

r 

r 

627 

-12.8 ±3.5 

-12.4 ± 3.5 

6.5 ± 10.9 

8.0 ± 0.4 

-1.5 ± 10.9 

3.1 

-0.5 ±3.5 

CfA4i 

r 

r 

1035 

-9.1 ± 3.4 

-4.9 ±3.3 

25.1 ± 9.4 

8.0 ± 0.4 

17.2 ±9.4 

3.1 

5.6 ±3.1 

CfA42 

r 

r 

1035 

-12.5 ±3.4 

-9.0 ±3.4 

13.5 ±9.3 

7.9 ±0.4 

5.6 ±9.3 

3.1 

1.8 ±3.0 

CSPl 

r 

r 

309 

-7.8 ±3.5 

-6.4 ±3.6 

1.7 ± 10.2 

9.6 ±0.2 

-7.9 ± 10.2 

3.2 

-2.5 ±3.2 

SNLSO 

r 

r 

200 

-8.2 ±3.5 

-7.2 ± 3.5 

23.3 ±3.2 

20.6 ±0.6 

2.6 ±3.3 

3.0 

0.9 ± 1.1 

SNLSl 

r 

r 

272 

-5.3 ±5.2 

-4.7 ± 5.2 

23.5 ±5.7 

20.6 ±0.6 

2.8 ±5.7 

3.0 

0.9 ± 1.9 

SNLSl 

r 

r 

264 

4.0 ± 5.3 

4.0 ±5.3 

25.2 ±6.1 

20.6 ±0.6 

4.6 ±6.2 

3.0 

1.5 ± 2.1 

SNLS3 

r 

r 

569 

-3.9 ±3.4 

-4.3 ±3.4 

20.7 ± 5.2 

20.6 ±0.6 

0.1 ± 5.2 

3.0 

0.0 ± 1.8 

SDSSO 

r 

r 

3300 

-11.0 ±3.2 

-12.5 ± 3.2 

-4.8 ±0.9 

-2.9 ±0.7 

-2.0± 1.1 

3.3 

-0.6 ±0.3 

SDSSl 

r 

r 

1525 

-9.5 ±3.2 

-10.3 ± 3.2 

-5.4 ± 1.3 

-2.9 ±0.7 

-2.5 ± 1.5 

3.3 

-0.8 ±0.5 

CfAl 

i 

I 

19 


-4.7 ± 22.3 

138.2 ±67.0 

87.6 ±3.8 

50.6 ±67.1 

1.9 

26.1 ± 34.7 

CfA2 

i 

I 

117 

2.6 ±4.9 

3.0 ±4.8 

64.4 ± 25.8 

87.6 ±3.8 

-23.2 ± 26.1 

1.9 

-12.0 ± 13.5 

CfAS 

i 

I 

319 

9.2 ±4.4 

7.9 ±4.4 

88.4 ± 18.6 

78.4 ± 6.0 

9.9 ± 19.6 

1.9 

5.1 ± 10.1 

CfAK 

i 

I 

557 

-1.5 ±3.6 

-0.8 ±3.6 

3.6 ± 12.9 

12.2 ±0.7 

-8.5 ± 12.9 

1.9 

-4.6 ±6.9 

CfA4i 

i 

I 

969 

-0.1 ± 3.6 

0.8 ±3.6 

17.4 ± 10.3 

12.2 ±0.7 

5.2 ± 10.3 

1.9 

2.8 ± 5.5 

CfA42 

i 

I 

969 

0.4 ± 3.6 

1.4 ±3.6 

17.4 ± 10.3 

12.2 ±0.7 

5.2 ± 10.3 

1.9 

2.8 ± 5.5 

CSPl 

i 

I 

306 

11.7 ± 5.4 

11.3 ±5.4 

11.4 ± 10.8 

14.4 ± 0.9 

-2.9 ± 10.8 

1.8 

-1.6 ±6.0 

SNLSO 

i 

i 

200 

-8.1 ± 5.2 

-6.6 ±5.3 

22.6 ±5.1 

16.2 ± 1.0 

6.3 ±5.2 

1.9 

3.3 ± 2.7 

SNLSl 

i 

i 

272 

1.4 ± 5.3 

1.8 ±5.2 

23.8 ±6.7 

16.2 ± 1.0 

7.6 ±6.8 

1.9 

4.0 ± 3.6 

SNLSl 

i 

i 

264 

3.4 ± 5.3 

3.2 ±5.3 

31.4 ±5.7 

16.2 ± 1.0 

15.2 ± 5.8 

1.9 

7.9 ±3.0 

SNLS3 

i 

i 

569 

-0.8 ±3.4 

-1.7 ± 3.4 

24.6 ± 4.4 

16.2 ± 1.0 

8.4 ±4.5 

1.9 

4.4 ± 2.4 

SDSSO 

i 

i 

3300 

-5.1 ± 3.4 

-5.6 ±3.4 

-2.9 ± 1.0 

-7.2 ±0.4 

4.3 ± 1.1 

1.8 

2.4 ± 0.6 

SDSSl 

i 

i 

1525 

-8.6 ±3.4 

-9.7 ± 3.4 

-4.6 ± 1.4 

-7.2 ±0.4 

2.7 ± 1.4 

1.8 

1.5 ±0.8 

CfAl 

z 

I 

17 


2.4 ± 20.7 

-22.3 ±54.8 

-105.1 ± 6.0 

82.8 ±55.1 

1.9 

42.7 ± 28.4 

CfA2 

z 

I 

116 

7.7 ± 6.6 

8.2 ±6.6 

-127.2 ± 21.3 

-105.1 ± 6.0 

-22.1 ± 22.1 

1.9 

-11.4± 11.4 

CfAS 

z 

I 

363 

2.1 ± 5.5 

8.9 ±5.5 

-80.7 ± 17.3 

-101.7 ± 10.3 

21.1 ± 20.1 

1.9 

10.9 ± 10.4 

CfAK 

z 

I 

675 

-5.6 ± 7.0 

0.1 ± 7.0 

-189.7 ± 10.5 

-168.0 ± 15.5 

-21.7 ± 18.7 

1.9 

-11.7± 10.0 

CfA4i 

z 

I 

1117 

-8.4 ± 6.8 

-0.5 ± 6.7 

-168.7 ± 9.4 

-168.0 ± 15.5 

-0.7 ± 18.1 

1.9 

-0.4 ± 9.7 

CfA42 

z 

I 

1117 

-7.9 ±6.8 

0.0 ±6.7 

-168.7 ± 9.4 

-168.0 ± 15.5 

-0.7 ± 18.1 

1.9 

-0.4 ± 9.7 

CSPl 

z 

I 

337 

4.9 ± 7.2 

9.2 ± 7.2 

-175.3 ± 11.2 

-165.8 ± 15.4 

-9.6 ± 19.0 

1.8 

-5.3 ± 10.6 

SNLSO 

z 

z 

158 

7.7 ± 4.1 

11.6 ±4.1 

30.2 ±4.9 

17.4 ± 2.0 

12.8 ± 5.3 

1.6 

8.0 ± 3.3 

SNLSl 

z 

z 

210 

5.6 ± 5.7 

7.3 ±5.7 

30.4 ±8.2 

17.4 ± 2.0 

13.0 ± 8.5 

1.6 

8.1 ± 5.3 

SNLSl 

z 

z 

264 

10.8 ±5.7 

11.4 ±5.7 

40.0 ±5.8 

17.4 ± 2.0 

22.6 ±6.1 

1.6 

14.1 ± 3.8 

SNLS3 

z 

z 

468 

4.3 ±4.2 

6.2 ±4.2 

26.6 ±5.4 

17.4 ± 2.0 

9.1 ± 5.7 

1.6 

5.7 ± 3.6 

SDSSO 

z 

z 

3295 

16.4 ± 3.9 

16.6 ±3.9 

37.3 ± 1.6 

27.5 ±2.9 

9.8 ±3.3 

1.6 

6.3 ± 2.2 

SDSSl 

z 

z 

1525 

16.5 ±3.9 

16.8 ±4.0 

35.5 ±2.1 

27.5 ±2.9 

8.0 ±3.6 

1.6 

5.2 ± 2.3 


Offset and slope differences between each system and PSl. The second and third column show the PSl filter and comparison system 
filter. The fourth columns shows the total number of overlapping stars between 0.3 < g — i < 1.0. The fourth column shows the offset 
found when using the HST Calspec library to find a nominal calibration offset between the systems. These offsets should be added each 
survey’s magnitudes to agree with PSl. Columns 5-10 show both the calibration offset and difference in predicted slopes when using the 
NGSL library for the spectral transformation. Columns 9 and 10 show how to convert the difference in predicted and recovered slopes of 
the transformation to a change in the mean wavelength of the comparison filter. The offset given when using HST Calspec library is for a 
color range of 0.35 < g — i < 0.55 while the offset and slope when using using the NGSL library is for the color range of 0.3 < g — i < 1.0. 
CfAl values are missing due to lack of comparison stars for a small color range. 
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are expected to be around 3 mmag in each filter. Since 
most of the comparisons between PSl and other surveys 
involve stars from multiple fields, we estimate the sys¬ 
tematic uncertainty in the PSl Ubercal photometry to 
be the quadrature sum of an error floor of 3 mmag in 
addition to the filter-specific uncertainty (5,4,4,4 mmag) 
divided by the square root of the number of fields. 

The systematic uncertainty associated with the spec¬ 
tral library used for the synthetic transformations can 
be determined by comparing the transformations using 
multiple independent spectral libraries. We analyze stel¬ 
lar spectra from 5 different libraries: the HST C alspec 
library (version 005]^ the NGSL spectral library (Heap 


& Lindler|2007), the PNGS s pectral librar y (Pickles et al. 
in prep), the ‘Pi ckles Atlas’ (|Pickles|199 8|) and the Gunn- 
Stryker library (Gunn & Strykerj 1983). We must make 
specific cuts for each library to properly alleviate system¬ 
atic biases and optimize the most consistent comparisons. 
For the HST Galspec library, the NGSL library and the 
INGS library, we explicitly only include solar analog stars 
(—0.6 < [0/H] < 0.3). This cut allows for the most di¬ 
rect comparison between the majority of main sequence 
stars observed and the stars from the synthetic libraries. 
For the HST Calspec library, we only include standards 
that have been recently calibrated with WFC3 so that 
only the latest updates are used. Because of known flux 
biases with distance from slit center in the NGSL obser¬ 
vation^^ we only include standards that were observed 
within 0.5 pixels of the slit center. There is significant 
overlap between the INGS and NGSL libraries as the 
INGS library uses many NGSL spectra and includes its 
own correction for slit-loss in the NGSL observations. 
We do not exclude any spectra from the Pickles library 
or the Gunn-Stryker library. 

Examples of differences between these libraries are 
shown in Fig. where we find the synthetic transforma¬ 
tions between PSl and SDSS as well as PSl and CfAS. 
We only show the transformations for a color range of 
0.3 < gpi — ipi < 1.5 as that is the only part of the main 
sequence where many stars can be found and there is the 
smallest amount of scatter in the transformation. We see 
that there is clearly less of a dependence on the spectral 
library for the transformation between PSl and SDSS as 
the mean effective wavelengths and wavelength range of 
these filters are much closer. However, the dependence 
can be quite significant when comparing filters between 
PSl and GfAS; separation between the mean effective 
wavelength of filters from these systems can be > 200 A. 

Since we chose to use the HST Galspec library as our 
primary spectral library, we determine the systematic 
uncertainty in the synthetic transformation between sys¬ 
tems by instead using the HST NGSL library. As the 
NGSL library was also acquired using HST, it is a useful 
comparison to the HST Galspec library. Included in the 
uncertainties given in Fig. 3, we find discrepancies in the 
offsets when using these two libraries to be in the range 
up to 4 mmag. Differences between the HST Galspec 
library and other libraries, like the Gunn-Stryker library 
or the Pickles library, can be significantly larger. We are 


Differences in absolute flux of Calspec standards between ver¬ 
sion 003 and version 005 are on the 8 mmag level 

More information found here: 

https: / / archive.stsci.edu / prepds/stisngsl / 


limited to use a small color range when comparing sys¬ 
tems due to the small amount of standards in the HST 
Galspec library. We can increase this color range by using 
the NGSL library to extend from 0.35 < gpi — ipi < 1.0. 
However, since we do not have multiple HST Galspec 
standards in this larger color range, we cannot use two 
libraries to assess what an appropriate systematic uncer¬ 
tainty would be for this range. 

The impact of dust depends on whether there is a dif¬ 
ference in the dependence on color between the dust red¬ 
dening vector and temperature vector. This issue can be 
largely removed by accounting for the extinction in the 
field that each star is located in. The uncert ainty in the 
extinction values is discussed in S14. From [Schlafly &[ 


Finkbeiner (2011), we can safely assume that since the 
stars used m this analysis are removed from the galactic 
plane, all of the dust represented by the MW extinction 
value is between us and the stars used for this analy¬ 
sis, and any uncertainty in this assumption is included 
in the total systematic uncertainty of the extinction val¬ 
ues. We propagate the uncertainty of extinction values 
through the Supercal process and find that for an extinc¬ 
tion oi E{B — V) = Q.l mag, we can expect systematic 
uncertainties of 2, 0.5,0.5,1.0 mmag in griz respectively. 


TABLE 5 

Main Systematic Uncertainties in the Supercal Process 


Systematic: 
Spectral Library 
Dust 

PSl Ubercal 
PSl Non-linearity 


Uncertainty 
0 < ir < 4 mmag 
0.5 < X <2 mmag 
3 < ir < 6 mmag 
0 < ir < 3 mmag 


Top systematic uncertainties in the Supercal process. A range of 
uncertainties is given here, while the individual systematic 
uncertainty for each system/bandpass is included in the 
discrepancies shown in Figure 3 and given in Table 4. 


We investigate in Fig. 5 the systematic uncertainties 
due to non-linearity of the PSl stellar magnitudes. Here 
we remove the nominal j3 offset shown in Fig. 3 for each 
survey to probe possible biases with magnitude. Over¬ 
all, we choose to quantify the possible non-linearity due 
to PSl by finding the trend for each filter when com¬ 
paring with SDSS. This is roughly 0,2,3,5 mmag for 
gpirpifpi^pi respectively over a 3 mag range from ~15- 
18 mag of each filter. It is unclear how much of this 
non-linearity to attribute to PSl, since comparisons with 
SNLS and other surveys do not appear to favor the same 
biases with magnitude that SDSS does. Therefore we 
conservatively claim that half of this non-linearity is due 
to SDSS and half to PSl. From this plot, we also see 
significant trends at the faint magnitudes for the low- 
z surveys. Likely, these trends are indicative of either 
selection biases or poor photometry of fainter stars in 
the low-z surveys. These biases may also affect the SN 
photometry, though this will require analysis in a future 
study. 

We limit the systematic uncertainties to ones that are 
based on PSl and the methodology used here, and not 
of any of the other individual surveys. Exploring the 
comparisons between PSl and other surveys, many in¬ 
teresting pathologies can be found. For example, there 
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Fig. 4.— Left: Residuals of the synthetic transformations between PSl and SDSS when using different spectral libraries. The gpi — fpi 
color is used as the transformation color for these systems. The transformation with the NGSL library is used for the baseline trend 
and is removed for each case. The HST Calspec standards are represented by points on the plots. Right: Residuals of the synthetic 
transformations between PSl and CfAS when using different spectral libraries. The colors ScfAS ~ *P1 9pi ~ *P1 ^re used as the 
transformation colors for these systems. 


are certain fields of stars in the low-z surveys that appear 
to have a significantly different offset when comparing to 
PSl than the majority of fields of these surveys. One 
example can be seen in the comparison between i^pi and 
ffcsp in Fig- 6. While we find a negligible dependence of 
the relative offsets based on magnitude, we do see some 
dependence on MW dust reddening. However, this trend 
depends strongly on a handful of fields with high red¬ 
dening, and it is unclear if the trend is caused by the 
reddening or if those particular fields happen to be dis¬ 
crepant for other reasons. Overall we see on the order of 
5-10 mmag trends with brightness, dust or R.A./Dec. 
for CSP as well as many of the other low-z surveys. The 
5-10 mmag trends represents a rough estimate of the 
systematic uncertainties in these samples. 


4.3. Comparing SN Photometry 

A more direct method to determine calibration dis¬ 
crepancies between various filters with the ultimate goal 
of consistent SN photometry would be to compare the 
photometry of SN observed by multiple surveys. Unfor¬ 
tunately, this approach is limited by both the statistics 
and the methodology — comparing photometry of SN is 
much more difficult as the SNIa features are relatively 
deep and broad relative to stellar a bsorption, yet narrow 
compared to the width of a filter. Mosher et al.| (20121 
compared photometry of 9 SNIa observed by both SUSS 
and CSP, though argued that 4 of these SN could not be 
used as their cadences were not satisfactory or the SN 


were not typical SN la, which limited the comparison. 
Of the remaining 5 SN, Mosher et al. (2012) found that 
the photometry of each individual SN agreed to better 
than 1%, though the sample scatter is up to 8 % and 
there i s no s ingle, consistent offset for all SN. |Betoule| 
et al. ([201^ performs a similar comparison to that of 


iVlosher et al. (2012) when trying to compare CfA3 and 
CSP. Their sample has 17 SN, and they find offsets in 
BVri of 5 ± 4,-9 ± 3,24 ± 4,3 ± 12 mmag with scatter 
of 42, 21, 39,49 mmag respectively. No systematic uncer¬ 
tainties are given for this approach, and we find using Su- 
percal that these differences are 10±8, 3±9,4±9, —4±6. 
A larger sample is needed to understand if there are er¬ 
rors beyond those seen in calibration, such as those due 
to image subtraction or astrometric registration. An al¬ 
ternative way to assess whether Supercal improves the 
consistency of these two surveys is to look at the agree¬ 
ment in the measured distances for the same SN using 
photometry from separate surveys. We find without Su¬ 
percal, the difference in distances between the CSP and 
CfAS surveys for 23 SN is 0.031 ± 0.026 mag and af¬ 
ter Supercal the difference is —0.012 ± 0.026 mag. More 
statistics will be needed for a more robust diagnostic. 


4.4. Effects on Recovered Cosmology 

We may choose to remove the calibration discrepan¬ 
cies found between all the systems so as to create a more 
uniformly calibrated sample. We are able to force all 
systems to be relatively calibrated in a consistent man- 
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Fig. 5.— Relative offsets between stars in the PSl sam¬ 
ple and other samples as a function of apparent brightness in 
ffpi^pifpi^ipi• The offsets shown are relative to the nominal Su- 
percal correction given in Fig. 3. Comparisons for each survey are 
shown, and in black the offsets for SDSS are highlighted to quan¬ 
tify potential systematic uncertainties from the PSl non-linearity. 
Vertical lines show where the PSl sample is cut on the bright end. 


ner with a given system or the average calibration from 
multiple systems. For the average solution, we use the 
calibration from PSl, SNLS and SDSS, which are the 
most recently calibrated systems, to determine a joint 
solution of the baseline-Supercal calibration. We use 
the offsets shown in Fig. 3, and weight them by the 
systematic uncertainties given in Table 1. After doing 
so, we find average offsets from the PSl calibration of 
—5 ± 3, —7 ± 3, —2 ± 2,6 ± 4 mmag where the uncer¬ 
tainties given are from both the systematic uncertainties 
given in Table 1 combined with the uncertainties from 
the Supercal process. 

After making the systems consistent, we can then re¬ 
determine the cosmological parameters derived from the 
full sample. The light curves of all SN are fit with the 
SNANA packag e jKessler et al . |200^ so that the fits are 
consistent with Betoule et al. pOldl); the same SALT2 
model and host-mass Hubble residual step of 0.06 mag 
are applied. Further information about light-curve fit¬ 
ting, distance bias corrections and host masses will be 
discussed in the next PSl cosmology analysis (Scolnic et 
al. in prep). The number o f SN for each system , af¬ 
ter quality cuts explained in iBetoule et al.l (20141), is: 
CfAl (10), CfA2 (18), CfAS (33), CfAK(58)7TTA4_l 
(34), CfA4_2 (9), CSP(32), SDSS (359), PSl (111), SNLS 
(235). Differences in the recovered cosmology, when in¬ 


cluding only CM B constraints from [Planck CollaboraT 


tion et al. (2014), are shown in Table 6 for the unaltered 


Fig. 6. — Relative offsets between stars in the CSP sample and 
PSl as a function of apparent brightness in ppi, Milky Way red¬ 
dening, R.A. and Dec. The offsets shown here reflect the values 
after the Supercal correction from Fig. 3 is applied for the gcsp 
filter. Points in red are excluded from the offset calculation as part 
of the iterative 3-a clipping procedure. 


calibration, as well as when the calibration is forced to 
agree with the SDSS, SNLS, PSl or the average calibra¬ 
tion. The weighting of each sample is determined by the 
error shown in Figure 3 for each system, and we also in¬ 
clude calibration errors from the SALT2 model and from 
the HST AB system as part of the weighting covariance 
matrix. We find that values for w may change up to 
Aw ~ —0.040 depending on which primary calibration is 
used. For the average Supercal solution based on SDSS, 
SNLS and PSl, the change in w is —0.026. The pri¬ 
mary cause of these changes is due to the uncertainty in 
the B and V bands which affects the low-z SN. For the 
average Supercal calibration, we show the mean Hubble 
residual for each SN subsample in Fig. 7 (top) and the 
difference in distances as a function of redshift when the 
Supercal correction is applied relative to when it isn’t 
in Fig. 7 (bottom). The Supercal solution appears to 
significantly improve the agreement between the low-z 
samples and the ACDM model. We find mostly positive 
Hubble residuals from the ACDM model due to the host 
mass correction and the distance-bias correction (for lat¬ 


ter correction, see Figure 5 in Betoule et al. (2014|); more 
discussion in upcoming Scolnic et al. in prep). 'I’he net 
change in distances for each subsample is up to 0.055 
mag for the low-z systems, but only 1% for the higher-z 
systems. 

The Supercal process has a small effect on the distance 
scatter of the joint sample. For the entire joint sample 
shown in Fig. 7, for each of the five cases shown in Table 
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6 , we find a relative of [0.0, —2.3, -1-6.7, —9.9, -1-0.3] 
respectively for 906 SN. However, if we only measure 
the improvement in for the SNe with 2 < 0.1, we 
see a more significant improvement in each case: 
of [0.0, -3.2, -10.7, -7.2, -11.2] respectively for 225 SN. 
The effect of Supercal is obviously much larger for the 
low-z sample as offsets between PSl, SDSS and SNLS 
are small and the calibrat ion of SDSS a n d SN LS has 
already been connected in Betoule et al. (20131. In a 


full cosmology analysis, SAL'1'2 should be retrained with 
the optimal Supercal solution and simulated b iases of 
distance resi duals with color should be removed (Scolnic 
et al. 2014b). These improvements should better show 
the iihpact of the Supercal corrections. 


TABLE 6 

Differences between the various surveys 


Primary Calibration 

Aw 

AO,m 

No-Correction 

0.000 ± 0.052 

0.000 ±0.017 

Supercal-Avg 

-0.026 ±0.051 

-0.005 ±0.015 

Supercal-PSl 

-0.040 ±0.055 

-0.015 ±0.015 

Supercal-SDSS 

±0.010 ±0.052 

±0.005 ±0.017 

Supercal—SNLS 

-0.032 ±0.050 

-0.007 ± 0.016 


Differences in the recovered cosmology when calibration is forced 
to agree with that of a particular survey. The no-correction value 
is given when no change to any calibration is made. Errors given 
represent the total systematic error of the best fit cosmology 
and include the uncertainties from the Supercal process for each 
filter/system. 
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5. DISCUSSION 

The accuracy of the Supercal technique presented in 
this paper depends on the number of stars common to 
PSl and the comparison system, as well as the differ¬ 
ence in effective wavelengths for the filters that are being 
compared. While PSl has covered nearly the limit of its 
observable sky, it would still be possible to increase the 
PSl statistics of bright stars by observing at shorter ex¬ 
posure times. This could help the absolute calibration of 
PSl as well as the Supercal process. It is no longer possi¬ 
ble for most of the other surveys analyzed here to obtain 
more observations of stars for other surveys, though some 
could work to release more of the data already acquired. 
A clear path to improved accuracy would be to redefine 
the absolute calibration of the low-2: surv eys on multi¬ 
ple standards observed by Landolt (1992), rather than 
BD-f 17. In the current scheme, the accuracy of the cal¬ 
ibration of the low-z survey is reliant on the accuracy 
of the observations of BD-l-17 by both Landolt (1992) 
and RST. As BD-l-17 is a binary star, this strategy is 
sub-optimal. 

Furthermore, a complete understanding of the impact 
of these relative differences between surveys can be un¬ 
derstood only after the light-curve models are retrained 
with the new calibration. As the light-curve training is 
primarily based on low-z SN, any systematics in the cali¬ 
bration of these SN has an appreciable effect on the fitted 
SN distances. 

The Supercal process is primarily used here to deter¬ 
mine zeropoint offsets between the surveys, and does not 
correct for the color terms as given in Table 4. This can 
be done if there were more RST standards that were red 
{g — i > 1.0). Red standards are particularly helpful for 


Fig. 7. — (Top) Mean relative Hubble residual offsets to the 
ACDM model for SN observed by each system analyzed. The 
points in black represent the residuals when no change to the 
the calibration systems are made, and in red, they represent the 
offsets after we apply the average Supercal correction. (Bot¬ 
tom) The change in distances due to forcing all calibration to 
agree to the average Supercal solution. The values are such that 
A/r = psupercal — MNoCorr.- SN from different samples are repre¬ 
sented by different symbols and colors as indicated in the legend. 

the PSl, SDSS and SNLS surveys as SN at z > 0.5 have 
g — i > 1.0. 

Instead of waiting for more observations of standard 
stars to be obtained by HST, we can apply the same 
Supercal process discussed above to catalogs of stars ob¬ 
served by HST. Photometry from HST is ideal to use for 
cross-calibration as it is most closely tied to the HST Cal- 
spec AB system, it is photometrically stable and closely 
monitored. The difficulty is that observations taken with 
HST typically are done as part of small programs that 
have very different objectives, saturation limits and se¬ 
lection effects. In order to bypass the various differences 
in the small HST programs, we co uld focus on one very 
large calibrated field: COSMOS (|Scoville et al.||2007|). 
Only imaging in a single band, ACS/F8I4W (roughly f), 
was obtained for the entire field and while the relative 
calibration across the field is better than 3mmag, there 
was no effort to define the absolute flux of the observa¬ 
tions. This is not the case for other observations taken 
of the COSMOS field, so we can use other observations 
acquired from the Hubble Legacy Archive to determine 
the HST zeropoint of the entire field. The expected ze¬ 
ropoint does not differ across the sky and any change of 
the zeropoint with time is negligible. We find a zeropoint 
when calibrating to the AB system of 0.268 ± 0.003 mag 
such that F814Wab = F814Wcosmos — 0.268 mag. Ap- 
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plying the Supercal process, we find the calibration offset 
relative to PSl of 9 ± 10 mmag with the i band compar¬ 
ison and 6 ± 14 mmag in the 2 band comparison. The 
errors are relatively large because the F814W filter has a 
relatively different wavelength range and mean effective 
wavelength than that of either i-pi or z^i. More work 
should be done to tie HST observations to PSl so it can 
be included in the Supercal process. 

A final question remains about what is the best way 
to use the Supercal corrections in future cosmological 
analyses which use SNla. A simplistic approach could 
be not to correct for any of the Supercal offsets, but in¬ 
stead to use them as a systematic uncertainty which one 
could propagate into the full systematic covariance ma¬ 
trix. The inherent assumption of this decision is that 
the current calibration of each survey is optimal, given 
some uncertainty. This assumption certainly appears to 
be false. Instead, we prefer to correct for the Super¬ 
cal offsets based on a weighted calibration of the mul¬ 
tiple surveys, as done for the average case in Table 4. 
This approach creates a homogeneous photometric cal¬ 
ibration between the different systems and best ties all 
the samples to the AB system. For samples that cannot 
be included in the Supercal analysis, due to insufficient 
numbers of stars or due to the lack of a public release of 
local standards, we think these samples should no longer 
be included in a cosmological analysis. 


6 . CONCLUSIONS 

In this paper, we have presented a new method for 
SN calibration. We have determined the relative SN- 
zeropoint offsets between different filters of different sys¬ 
tems and show their impact on the measurement of cos¬ 
mological parameters. We find that there may be system¬ 
atic discrepancies between the zeropoints of 1 — 2% which 
propagate to up to 5% systematic errors in SN distance. 
This systematic error in distance can then result in an 
average offset of 2 — 3% in w. The primary discrepancies 
found are in the B filters of the low-z systems. More 
work must be done to better understand these offsets. 
The systematic uncertainty of this approach are shown 
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